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1 Introduction 



The variational approximation method has been proved to be a convenient non- 
perturbative approximation method for both the quantum mechanical and quan- 
tum field theoretic systems. The method, which was originally developed for a 
non-perturbative evaluation of quantum mechanical energy eigenvalues, has been 
extended to an action principle by Jackiw and Kermaic, which enables one to com- 
pute the effective actions for time-dependent systems and also for quantum field the- 
ories. The usefulness of the method for quantum field theories, however, is limited 
by the fact that only the Gaussian functional and it's obvious generalizations can be 
used as trial wave functionals, which is the reason why it is often called the Gaussian 
approximation. Attempts to improve the Gaussian approximation by using more 
general non-Gaussian trial wave functionals have not been asJruitful as expected, 
due to the difficulties in solving the variational gap equationta. To overcome such 
limitations the authors of Ref. [3] has developed the so-called post-Gaussian expan- 
sion method, which enables one to compute the order-by-order correction terms to 
the Gaussian effective potentials for the time-independent systems. 

It is the purpose of this note to give a brief review of the variational approxi- 
mation method in quantum field theory, and to present a new systematic method 
of perturbative expansion around the Gaussian effective action based on the back- 
ground field method, which is also applicable to the time-dependent systems. 

The variational method in quantum mechanics is based on the fact that the 
ground state energy is the minimum of the expectation values of the Hamiltonian 
of the system. This can be written as a variational principle, 

6^ <^\H{Q,P)- >=0, (1) 

where H is the Hamiltonian of the system and > is an arbitrary state. If one 
starts from a trial wave function with variational parameters, then Eq.(|l|) provides 
a variational approximation of the energy eigenstates. 

This simple variational method can be generalized to include the information 
on the time-evolution of the system, by choosing the state to satisfy the constraints, 

< ^-l^- > = 1 
<^\Q\^> = q{t) (2) 
<*|P|^'> = p{t). 

Minimizing the expectation value of the Hamiltonian with respect to such states, 

(5* < >=0, (3) 

we obtain the effective Hamiltonian, 

H,ffiq,p)^<^o\HiQ,P)\^o>, (4) 

where |4'o > is the solution of Eq.(^. This effective Hamiltonian provides many 
useful informations on the quantum system. If one minimizes Hef / with respect to 
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the variations of q{t) and p{t), one obtains the ground state energy eigenvalue of 
the system. If one solves the Hamilton's equation, on the other hand, 



m = ^ (5) 

*) - 

the solution q and p give the quantum mechanical informations on < > 
and < \l/o|-P|^o >: respectively. 

If we set (7 = = p in Eq. (||) , we obtain the effective potential, 

Veffiq) = H^ff{q,p)\p=o- (6) 

If one starts from a trial wave function which satisfies Eq.(||), then one obtains 
the variational approximation of the effective Hamiltonian or the effective potential. 

This simple argument is in fact a part of the quanium action principle in the 
Schrodinger picture called the action principle of DiracrJa. This formulation can be 
generalized to the case of quantum field theory, which is the functional Schrodinger 
picture formulation of quantum field theorjo. 

2 Functional Schrodinger picture representation of quantum field the- 
ory and the variational method. 

2.1 Schrodinger picture representation of scalar field theory 
The scalar field theory is defined by the Lagrangian density, 

£ = i(^0-(/)7V)-H[(/.,7r], (7) 

and the equal-time commutation relations 

[4>{x),Tr{x')]^o=^,o = i6{x~x') (8) 
[0(x),0(a;')]^o=^,o = = [7r(a;),7r(a;')]:cO=x'o, 

where Ti,[(f>,TT] = ^vr^ + ^(0) is the Hamiltonian density. 

As in the quantum mechanical case, we choose the basis vectors for the Schrodinger 
representation to be the eigenstates of the field operator (j){x), 

(l3\ip{x) >= ipix)\ip{x) >, (9) 

which satisfy the orthogonality relation and the completeness relation, 

< ^1^' >= Sy{x) - ^'(.f)] (10) 
Vifilif X (p\ ^ I. 
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Then the state of the physical system is represented by the component, ^'[(p, t] = 
< (p\'^ >, of the expansion, 

>= J VLp\tp><tp\^ >, (11) 

and the operators are represented by the matrix elements, 

<(p|0|^'> - ip{x)S[v ~ f'] (12) 

which are obtained from the matrix elements of Eq. (^ . The time-evolution of the 
system is then determined by the functional Schrodinger equation, 

i^*[(^,t] = < (^|i/[0,7r]|v3' >< (^'l^- > 

= H[^,-t-^M^,t]. (13) 

This is a simple generalization of the Schrodinger picture representation of quantum 
mechanics to the case of the systems with infinite degrees of freedom. 

2.2 Free Scalar field theory 

To obtain physical information in the Schrodinger picture formulation one has to 
solve the functional Schrodinger equation ( [l3| ) . As an example of solving the equa- 
tion of motion, we consider the free scalar field theory described by the Hamiltonian 
density, 

H=i7r2 + i(V0)^ + imV. (14) 
By using the representation (|l^), the Hamiltonian of the system can be written as 

H=IJ - ^(^)(V^ - m^Mx)]. (15) 

Since the Hamiltonian has no explicit time-dependence, there exist stationary 
states satisfying the energy eigenvalue equation, 

H[^,-t^]np]^E'f[^], (16) 



where E is the energy eigenvalue. Since the Hamiltonian (|l^) is quadratic in ip{x), 
the Gaussian functional 

^[ip] — / d^xd^x'ip{x)G{x,x')(p{x') ^j^y-j 

is a solution of (p^). Substituting (|lj) into (^6|), we obtain. 

i[G(x,x')-¥'(^)G'(x,x')¥'(^') + ^(^)(-V'+m2)^(x)]4'[(p] (18) 



where the repeated arguments imply the integration over the variables. For this 
equation to be true for all eigenvalue ip{x), the following relations must be satisfied, 

E = ]-TrG =11 d^xG{x, x) (19) 
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G^{x,x') = i-V^ +m^)S{x~x'). 



Solving the second equation of ( |19[ ) via Fourier transformation, one finds the kernel 
of the Gaussian functional ( p^ ) is given by 

G-\x, x') = -i^ / d?pe^^<^-^')^^=, (20) 

which is the propagator of the free scalar field theory at equal-time. From this the 
eigenvalue E is found to be 
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E = 2 y J d^pV^T^, (21) 

which is the well-known zero-point energy of the free scalar field theory. 

We thus find that the ground state of the free scalar field theory is described 
by the Gaussian wave functional with the kernel given by the inverse propagator at 
equal-time. One can construct all the eigenstates of the Hamiltonian by considering 
the creation and annihilation operators. They are given by the Hermite polynomial 
functionals multiplied by the ground (vacuum) stateEl. 

For the interacting theory, however, we cannot solve the functional Schrodinger 
equation (^3|) in closed form. One thus has to use an approximation method to 
find solutions to the equation of motion. A natural way to evaluate an approximate 
solution to ( [T^ ) in the Schrodinger picture is to use the variational method, which 
is the subject of the next subsection. 

2.3 Gaussian approximation of (j)'^ theory 

To illustrate the variational approximation method in quantum field theory, we now 
consider the scalar (j)'^ theory described by the Hamiltonian density, 

w = ^-^ + i(w)^ + ^/^v+^0^ (22) 

and compute the Gaussian approximation of the effective potentiaEI. To do this we 
take as a trial wave functional the solution of a free field theory discussed in the 
last subsection, namely, the Gaussian wave functional, 

which satisfies the constraints, 

= 1 

< * I I * > = 0(f , t) (24) 

< * I TT I > = 7r(f,t). 
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Since the effective potential is the effective Hamiltonian with vanishing mo- 
menta, we set TT equal to zero and consider the space-time independent (j). Then the 
expectation value of the Hamiltonian with respect to the trial state(|^ is given by 

<H> = Jd^x (v{^) + h[^G"\x,x) ~ ]^^G{x,x')\s=s' 

+ V^^\^)G{xM + yl^(^)(0)(G(x,x))2) , (25) 

where Y(p, ^ i^4>^ + ^^'^ is the classical potential and = (i|)"^(<^)- 

The variational equation, 

Sg<H >= 0, (26) 
yields the variational gap equation, 

^G-^{x,x') = [-V^ + v'^^\$) + ^v'^^H4>)Gix,x)]S(x-x'). (27) 
Substituting(^^ into ( p5| ) one obtains the effective potential, 

VeffW = + ^4>' + ^G-\x,x) ~ ^G(f,f)G(x,x). (28) 



The first three terms Eq.(p8[) constitute the one-loop effective potential, and the 
last term contains the higher-loop contribution, which shows the non-perturbative 
nature of the approximation. 

Using the Fourier transformation method we can solve the gap equation(^7|): 

(2^)^ 

where the effective mass m is defined by 



(27r) + 



2 _ 2 , 72 , t /■ 



+ -4)' + -h / ^ , (30) 

2"^ 2 7 (27r)^VP^T^ 

Substituting these results into ( p8| ) yields the effective potential: 

- — /2 M H ^ In — , 31) 

4 ^ ^ 647r2 V M2 2/ ' ^ ' 

where 

= /i + -^1 + -<r m^hiM) + ^ In — ^, 32 

^ 2 2^ 2 ^ ^ 327r2 M^' ^ ' 

and, /i and /2(M) are given by 



The effective potential (31) has two divergences, Ii and l2{M), and thus we 
need to renormahze the theory. The renormaUzation conditions can be written as 
[6], 



dVr- 



eff 



dm? 

which, in our case, yields 



ITlQ 



/4 



d(rr?Y 



1 



/4 
1 



A ^ 2^ 
1 1 
A 



(34) 



(35) 



-MM)- 



Substituting (pa) into (pl^), we finally obtain the renormalized effective potential. 



4 ? ? 

/ 7\ ™ 72 Mb 2 



64^ I ''m2 



where 



2 , 



Afl , rn^ 

77 In — ^. 

327r2 M2 



(36) 



(37) 



This result is qualitatively equivalent to the large- iV approximation (where N is 
the number of the scalar field componentsP, and shows the non-perturbative nature 
of the approximation. The result ( ^6| ) can be used to study the quantum nature 
of the theory, such as the symmetry breaking phenomena^. The method can be 
generalized to the case of the field theories in curved spacra, to the finite tertmerature 
quantum field theorp, and also to the systems out of thermal equilibriuntJ. 



2.4 Fermionic theory 

Finding the Schrodinger picture representation of fermion fields is non-trivial since 
the equal-time anti-commutation relation, in Hcrmitian representation. 



{'lpa{x),llj0{x')}^o^^o' = 5a(3S{x-x' 



(38) 



implies that the field operators ip^s are their own momentum conjugates. In the 
non-Hermitian representation, on the other hand, the anti-commutation relation is 
given by 

{'ll;a{x),1pl{x')}.^0=^0' ^ 6al3S{x - x'), (39) 

but now the eigenstates of field operators do not form an orthonormal basis. One 
can show, however, that^there exists a one param£ter family of equivalent represen- 
tations of fermion fielda. Floreanini and Jacki'vciJ, in particular, have found that 
the representation. 



(40) 



se{x)' 
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where d{x) is a complex Grassmamr variable, satisfies the anti-commutation relation 
@. 

Using this representation one can follow the same procedure as in the case of 
the scalar field theory to evaluate the Gaussian approximation of the Fermionic field 
theories. 

This method has been applied, to the cases of (l+l)-dimensional and (2+1)- 
dimensional Gross-Neveu model£j, and found that the results are qualitatively 
equivalent to those of the large- approximation. It has also been applied to the 
(2-|-l)-dimensional Thirring model, and in this case the GaussiaiL result has been 
found to be qualitatively better than the large- iV approximatioio. The Gaussian 
effective potential shows the existence of the symmetry breaking in the theory, while 
the large- does not. The method has also been-Jound to be effective in determining 
the existence of the bound states in the theor\til. 



3 Beyond Gaussian approximation? 

We have shown that the Gaussian approximation based on the functional Schrodinger 
picture formulation is a simple generalization of the quantum mechanical varia- 
tional method, and provides a convenient non-perturbative approximation method 
for various quantum field theories including the systems out of thermal equilibrium. 
However, the method has a serious drawback in that, unlike the quantum mechan- 
ical cases, one can compute the expectation value of the quantum field theoretic 
Hamiltonian only for the Gaussian wave functional and its simple generalizations. 

If one could develop a systematic method to compute the correction terms to 
this approximation, therefore, the variational method would become a very powerful 
technique in obtaining physical information from quantum field theories. In fact 
there have been some attempts in this direction recently. 

One of such attempts, known ^as the post-Gaussian approximation method or 
the variational perturbation theorp, is based on the observation that the Gaussian 
wave functional 

^^^^Q-kv(S)G-\x-x')v{S') ^ ^ = <^_0 (41) 

is a vacuum state for a free theory with mass given by 

+ (42) 

where g{p) is determined by the variational gap equation ( |29| ) and (^0[). 
One can construct all the eigenstates of this free theory: 

*„[(^] =i/„[^]*o[(^] , n = 0,1,2,- •• (43) 

where Hn[<p] is the Hermite polynomial functional. One splits the Hamiltonian as 

H = Ho + Hi (44) 

where Hq is the free Hamiltonian with the mass m given by Eq.(^2[) and Hj is the 
remaining part of H, and then applies the well known perturbation theory with the 
states (^) taken to be the basis for the perturbation. 
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The effective potential is the expectation value of the Haniiltonian with the 
constraint 

< *o|0|^'o >= (45) 

where |\l/o > is the ground state of the system. 

Taking the eigenstates ( p3| ) to be the basis for perturbative expansion, the 
effective potential can be written as 

oo / 1 \" 

V,ff{^)^Eo{$) + J2<'^o\'^i[^ T^Hj) \^o>, (46) 

^0 \Eo-Ho J 

where the first term is the Gaussian effective potential and the second term repre- 
sents the perturbative correction tcrmfl. This method has been used to compute 
the first non-trivial contributipii to the Gaussian effective potential for the (0-1-1)- 
dimensional (p^ theory by Jeoitj and the (1+1)- and (2+l)-dimensional (j)^ theories 
by Cea and TedesctEl. 

This method provides a conceptually simple procedure to compute the order- 
by-order correction terms to the Gaussian effective potential for time-independent 
systems. For the time-dependent systems such as those out of thermal equilibriun£3, 
however, we need a method that can deal with the effective action, which contains 
information on the time-evolution of the system. It is the subject of the next 
section to introduce a new systematic, expansion method for effective action, which 
was developed by Geon Hyoung Lei£l. 



4 The Background field method 

We now introduce a systematic method which enables one to compute the order by 
order correction terms to the Gaussian approximation of the effective action, that 
contains the information on the time-evolution of the system. This new method is 
based on the background field method. So we will start to give a brief review on 
the background field methocLD. 



4-1 Background field method 

The generating functional for the connected Green's functions, is defined by 

the path integral. 



where S[(j)\ is the classical action, and the integral convention, Jcf) = J d"'xJ{x)4){x) 
in n-dimensional space-time, is used in the exponent. The effective action T[<j)\ is 
then defined by the Legendre transformation, 

TM = W[J] - Jip, (48) 

where (p{x) is the vacuum expectation value of the field operator, 

^{x) =< 0(x) >,7- ^^W- (49) 
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The background field method is one of the most effective ways of computing the 
effective action in the loop expansion, and is obtained by introducing the background 
field B{x) to the classical action: 

S[cl)]^ S[cb + B]. (50) 
This defines a new generating functional, 

^iW[J,B] = N j ']:,(l)(,^S[4,+B]+iJ,p 

= N j p^'e^sM+'-'I^'-s]. (51) 

Eg.(pTl) implies that this new generating functional is related to the old one by 

W[J, B] = W[J]-JB, (52) 

and the vacuum expectation value of the field operator in the presence of the back- 
ground field is given by, 

5W 

^= — = ^-5. (53) 

Now consider the case when the vacuum expectation value of the field operator 
in the presense of the background field, (p, vanishes: 

5W 

jj-^ = 0. (54) 

This determines the configuration of J which extremizes the generating functional 
W[J,B]: 

J = J[B]. (55) 



Then one can easily show from Eq. (|5^) that the extremum value of the generating 
functional W becomes the effective action: 

W[J, B]\^=o = W[J{B)] - J{B)B - r[B]. (56) 

The simpler way to compute the effective action T[B] is to consider the Legendre 
transformation of W[J, B]: 

t[ip,B]^W[J,B]- Jip, (57) 

which generates the one tp-particle irreducible diagrams in the presence of the back- 
ground field. If we set = in Eq.(|57|), we find that 

f[^,B] = W[JiB),B]-Jip 

= r[B]. (58) 

This implies that the effective action r[B] consists of one i^-particle irreducible 
bubble diagrams with no external ip-lines. Thus one can write the effective action 
as 



(p=0 
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(59) 



To illustrate how the method facilitates the computation of the effective action, 
we now consider the scalar 0* theory as an example. The Lagrangian density of (j)'^ 
theory can be written as 

m - -^Hd' + m')<^ - ^d^' ^ -l<i>{^-')4> - (60) 

where D is the free propagator defined by 

V = i{d^ + (61) 
The shifted Lagrangian is then given by 



(62) 



where K is the new propagator, 

K^^{^^ + ^J:' + ^B^r\ (63) 

which includes a part of the interaction effect through the term proportional to 
B^ . The generating functional in the presence of the background field can then be 
written as, 

(64) 

The effective action for the scalar i/)^ theory now becomes 



gjr[i3] ^ ^tW[J,B] 



= Ne^''^ ""-^^ det{-K)Uo[B], (65) 

</3=0 



where 



(66) 



The coefficient of /o[-B] in Eq.(|65|) represents the one- loop effective action, and 
Io\B] represents the higher loop contributions to the effective action. The first 
factor of Eq.(|6^) constitutes a shift operator, which docs not change the extremum 
value of W[J, B]. Thus this factor, which generates the tadpole diagrams, does not 
contribute to the effective action , and we can neglect it . One can therefore write 
h[B] as 



(67) 

which implies that the higher loop contributions to the effective action consist of 
one (^-irreducible bubble diagrams with no external (^-lines. One can compute the 
higher loop contributions by the Feynman rules: 
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-iXB 



= -i\ 



3-vertex 



4- vertex 



One can easily show that the two-loop contribution to the effective action con- 
sists of two bubble diagrams with no external line£Il: 




Comparing this with the conventional method of computing the effective action, 
one can easily convince oneself that the background field method provides a very 
simple and effective way to compute the effective action in the loop expansion. 



4-2 Expansion around the Gaussian effective action 

The reason why the procedure of computing the effective action in the background 
field method is simplified is that the propagator used in the Feynman rules, Eg. (|63|) , 
includes a part of the interaction effect through the background field B{x), and the 
generating functional has been rearranged in such a way that the zeroth order term 
in the expansion, the coefficient of IqIB] in Eq.(|65|), is the one- loop effective action. 

We now want to proceed one step further and rearrange the effective action 
in such a way that the propagator used in the expansion would include the non- 
perturbative effect and the zeroth order contribution would be the result of the 
non-perturbative approximation, such as the result of the Gaussian approximation. 

Note that the Gaussian approximation consists of the cactus type diagrams 
where an internal line comes out of a point in space-time and goes back to the 
same point, such as the first term of Eq.(|68|). So we need to extract such diagrams 
out of the Io[B] part of the effective action. To do this we examine the functional 
derivatives of the generating functional e^-^'^-^ with respect to J{x) : 

G^yJye^'y'^y''' = (GJ),e5^^-^ (69) 

((GJ)2 + G,,)e^^«-' (70) 



6 1 7 r 
g 2 ^y^v 

SJx 



5Ji 
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^e^^-^-^^ = ((GJ)^. + 6G..(GJ)^+3GL)e^'^^^ (72) 

which would appear in the expansion of /o[^]- Note that these derivatives, which 
generates 2-loop and higher order contributions to r[B], generate the cactus type 
diagrams as can be seen in Eqs.([7l|) and (|7^). 

To extract such diagrams out of Io[B], we define new functional derivatives so 
that the following relations satisfy: 

i^yei'^' = {GJ)lei'^' (73) 



Comparing Eqs.(|7l|) 
must be defined by 







(74) 


H) and Eqs.dxl)- 


djl), wc find that the 


primed derivatives 


^5Jl' ~ SJi 




(75) 






(76) 



Note that the primed derivatives operated on the functional e'^'^^'' do not generate 
the cactus type diagrams that contribute to the Gaussian effective action. This 
implies that by using the primed derivatives one may in fact rewrite Io[B\ in such 
a way that all the cactus type diagrams would be included in the zeroth order term 
in the expansion. 

To do this we need to find a new Green's function G such that 

where a normalization constant N' and a function A are to be determined. Using 
the definitions (75) and ([76|), one easily finds G, N' and A that satisfy Eq.(|77|): 

N' = det(-G)^det(-G-l + yG,,)5e'^«'^ (78) 
A = ^BGxx, (79) 

^xy ^ ^xy ' 2 ^xx^xy 

= T^xy - -i^BlSxy + -:^Gxx5xy (80) 

Using these results we can write the generating functional in the presence of the 
background field as 

e''^^-'^''^ = det(-G)ie^^^"^-'5^^'+'^«- 
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The effective action can thus be written as 



det(-G)^ei^^" 



(82) 



where 



I[B] 



exp 



XB 



,XB , 5^ iX , 5'^ ... 1 ,r 1 



(83) 



(p=0 



In Eq.(g^) we have used the fact that the shift operator, the tadpole part of the 
first line, does not affect the extremum value. Note that Eq.(pq), which determines 



the Green's function Gxy, is exactly the variational gap equation written in the 



Minkowski space notation, and therefore the coefhcient of I[B] in Eq.(|8^) gives rise 
to the Gaussian effective action as we have alluded to above. Since I[B] can be 
expanded as a power series in A as is done in the conventional background field 
method, we have the perturbative expansion of the effective action around the 
Gaussian approximation. 



The structure of I[B] is the same as that of Iq[B] in Eq.(67), except that the 
propagator is now replaced by Gxy and the functional derivatives are replaced by 
the primed derivatives. Thus the I[B] contribution to the effective action consists 
of one (p-particle irreducible bubble diagrams with no external (^-lines, and without 
the cactus type diagrams. The contributions from I[B] can therefore be evaluated 
by the Feynman rules: 
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We have thus estabhshed the systematic method of computing the order-by- 
order correction terms to the Gaussian effective action. 

One can easily show that the order contribution of /[B], for example, consists 
of the diagrams, 




Note that the first diagram of Eq. ( |68| ) is a cactus type diagram and does not appear 
in the expansion of I[B]. 



4-3 Effective potential 

If we consider the space-time independent background field in the above formulation, 
we obtain the effective potential defined by 

Since Gxy is a function of {x — y) for a constant _B, one can solve Eq.(|80|) by Fourier 
transformation, 

G{x-y)^ J ^gip)e^P(--y). (86) 
The Green's function Gxy in the momentum space can then be written as 



dip) = TT- — 2T' (87) 
«(p^ — rriQ) 



where 

Note that Eq. ( p7| ) is the variational gap equation of Gaussian approximation written 
in the Minkowski notation. 

Up to the order contribution of I[B], the effective potential is therefore given 

by 

Veff{B) = VG + Vp, (89) 

where 

Vg ^ ^B' + ^B' - |^(m2 - - ^B^ + *Tr log V^, (90) 

and 

Vp^^^B'G'~^^G\ (91) 



Vg of (|90|) is the Gaussian effective potential and Vp of (91) is the first non-trivial 
contribution from I[B], which is represented by the Feynman diagrams shown in 
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(gj). This rcsuh is the same as that of Cea and TedesccEl computed from the 
variational perturbation theory mentioned earlier. 

To illustrate the effectiveness of the perturbative correction method, we now 
consider the (0 + l)-dimensional case, which is the quantum mechanical anharmonic 
oscillator. It is convenient to write the Green's function in Euclidean space which 
is defined by 

p2 = -p|, dy = id''pE. (92) 



(93) 



Then G{x,x) in Eq. (|88|) can be written as 

dp 1 f dpE 1 11 



27ri(p2_TO2) J 27rp| + m2 \/2\/^' 
and the gap equation becomes 

= + ^B^ + -^=. (94) 
2 4Vm2 

The effective potential up to the order of I[B] becomes 

Veff{B) = VG + Vp, (95) 



where 

2 4! 2A^ "^2 



u^T A. 1,T o A„„ V 



and 



X2 \2 

Vp = -jj^B' (97) 

144to4 1536(V^)5 



We can compare this result with other approximation methods. For the com- 
parison of the effectiveness of various approximation methods it is convenient to 
compare the ground state energy eigenvalue of the anharmonic oscillator, which is 
defined by the minimum value of Vef / . The table below shows the errors of various 
approximation methods compared to the numerical results^. This shows that our 



Table 1: Errors of Eo{y? = 1) compared to the numerical result 



A 


1-loop 


2-loop 


Gaussian Approx. 


Improved Gaussian 


0.24 


1.4% 


0.05% 


0.006% 


0.0003% 


24 


38% 


56% 


1.1% 


0.37% 



method greatly improves the Gaussian approximation even at the first non-trivial 
correction level, for all values of the coupling constant. Note that, for strong cou- 
pling, the Gaussian and the improved Gaussian approximation display much better 
results than the loop expansion methods, which shows the non-perturbative nature 
of the approximation methods. 
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5 Conclusion 



We have shown that the functional Schrodinger picture formulation of quantum 
field theory provides a conceptually simple non-perturbative approximation method, 
namely the Gaussian approximation. It has been proved to be a reliable method to 
obtain non-perturbative informations on the symmetry breaking phenomena, exis- 
tence of bound states and time-evolution of the time-dependent systems with rela- 
tively little computational difficulties. The method, however, has a serious drawback 
in that it does not provide a systematic procedure to improve the approximation. 

Wc have shown that the background field method provides a solution to this 
limitation of the variational method. We were able to rearrange the effective action 
in such a way that the zeroth order contribution in the expansion of the effective 
action consists of the Gaussian effective action. This provides a new systematic 
expansion method which enables one to compute the order-by-order correction terms 
to the Gaussian effective action. 

The effective potentials for the (0 -I- 1), (1 -I- 1) and (2 -I- 1) dimensional (j)^ 
theory computed by the new method agree with those of the variational perturbation 
theory. It has been shown that the perturbative correction greatly improves the 
Gaussian approximation even at the first non-trivial correction level. Our method 
can also be applied to the time-dependent systems and can easily be extended to 
the cases of fermionic and gauge field theories. 

Although the effective action expanded around the Gaussian approximation for 
the lower-dimensional theories can easily be shown to be renormalizable, one has 
to study the renormalization procedure systematically for the (3 -|- 1) dimensional 
case. 
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